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Abstract 

We show that the dynamical behavior of the 3D Ising spin glass with Gaussian 
couplings is not compatible with a droplet dynamics. We show that this is implied 
from the data of reference that, when analyzed in an accurate way, give multiple 
evidences of this fact. Our study is based on the analysis of the overlap-overlap 
correlation function, at different values of the separation r and of the time t. 

In a very interesting paper Komori, Yoshino and Takayama discuss the dynamical 
behavior of the three dimensional {3D) Edwards Anderson (EA) Ising Spin Glass with 
Gaussian couplings. A number of very accurate numerical evidences are used to try to 
determine if the 3D EA model behaves in analogy to the Replica Symmetry Broken (RSB) 
solution of the Mean Field Sherrington Kirkpatrick model Q, or if its behavior is the one 
typical of coarsening of domain walls |J , in agreement with a picture of separated droplets 
I- 

The evidences presented in Q can seem at first view to have mixed signs. For sake of 
completeness let us try to summarize all the available evidence. 

As it is well explained in the paper |lj] some of the findings are obviously devastating 
for a droplet picture. The free energy barriers Bl found on a lattice of linear size L do not 
grow with L like L^, but with a logarithmic dependence B L oc log(L) (this was already 
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observed by Kisker, Rieger, Santen and Schreckenberg in 0). Secondly, at variance with 
the droplet picture, the width of the distribution of the free energy barriers instead of 
increasing with L turns out to be remarkably constant. These two features do not have 
any problem with the RSB picture. 

Two more findings are not in contradiction with any of the two exemplar behaviors. 
The relaxation process has an Arrhenius-type form, with the characteristic time of the 
form r = r exp (jf), and the decay of the energy density to its infinite time asymptotic 
value follows a power law. Both the droplet picture and the RSB scheme imply a power 
law decay of the energy density. 

The crucial point, as rightly pointed out in jj]], is the dynamical behavior of the overlap- 
overlap (q — q) correlation functions G(r, t) (computed starting from two different random 
spin configurations, considering overlaps at distance r and equal time t). A simple scaling 
of G as a function of |tW would indeed be a strong argument in favor of a droplet behavior, 
while evidence for a different scaling points to a RSB pattern. The analysis of JIJ claims to 
show a droplet like behavior of G(r,t). Here we show that, on the contrary, the numerical 
data exclude this simple scaling form. We show that they behave in the non-standard 
way already discussed in reference ||. We perform a careful data analysis, and show in 
detail where the analysis of |IJ] fails. In this way we show that the observed picture is fully 
compatible with an RSB behavior, and that none of the available numerical supports a 
droplet-like behavior: one more crucial piece of evidence, the one concerning the behavior 
of the overlap-overlap correlation functions, conjures with the observed behavior of the free 
energy barriers to exclude a dynamical droplet like behavior of 3D Ising spin glasses. 

We will try to make clear to the reader how delicate this analysis is, how even a careful 
treatment like the one of reference fl[] can misdirected by only apparent straight lines, and 
why we can get a completely safe evidence of a non-trivial scaling of G(r, t). 

As we have already said the crucial difference in the approach to equilibrium among the 
droplet approach and the RSB theory is in the behavior of the q — q correlation function 
during the approach to equilibrium. 

Let us define 



G(r) = lim G(r, t) . (1) 

t — >oo 

In the droplet model we have 

lim G(r) = q ± , (2) 

1 — >oo 

while in the RSB solution of the Mean Field theory (in the q ~ sector of the theory: we 
are analyzing this case since we start out from equilibrium on a very large lattice) we find [] 

G(r) « r - " , (3) 

and therefore 

1 A few details about our runs: we have analyzed four samples of a 3D Ising spin glass with Gaussian 
couplings, periodic boundary conditions, on a lattice of linear size L = 64. We have computed the overlap- 
overlap correlation function for runs of time length t n = 100-2™, where n runs from 1 to 13. We average 
the correlation function over the whole run: this gives a measured G that is numerically slightly different 
from the one of |jj but has the same scaling behavior. We have used random initial configurations. 
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Figure 1: The overlap-overlap correlation functions G(r, t) versus | for different r and t 
values. 



lim G(r) = . 



(4) 



In reference || we have studied this problem and we have concluded that \x ~ .5, and 
therefore that a droplet-like behavior is inconsistent with the results of numerical simula- 
tions. 

Reference Q argues that the numerical data for r > 4, T = 0.7, scale like 



G(r,t) = M 



(5) 



where = t 7 ^ , z(T = 0.7) ~ 8.71 and M is a scaling function. As we have already 
argued the validity of this analysis would be a good point in favor of a droplet like behavior 
of 3D Ising spin glasses (at variance with the many other evidences of that we have 
discussed before, that point against the validity of a droplet picture). Since this conclusion 
is different from the one we had reached in || we will analyze here its derivation in better 
detail. In the rest of the paper (in particular in figures [l], H and we use the value 
z(T=0.7) =8.71. 

The conclusion of [0] is based on the plot of the raw data for G(r, t) as a function of 
|£r for r > 4 (figure 5 of |TJ). So at first, mainly to exclude the possibility of programming 
errors, we compare our raw data (0, 0]) to the ones of figure 5 of 0. We plot in figure [1] our 
raw data for G(r,t) versus |: by comparing with [|TJ one can see they are not substantially 
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Figure 2: G(r = 4, t) versus I. These data are not on a straight line. 



different, so we can feel safe about the quality of the numerical data of both |l| and || 0. 
In figure [I] we have selected the same x and y scales than in figure 5 of JIJ , to make the 
comparison of the two sets of data easier. In order to make a visual comparison easier we 
have also selected time values as close as possible to the ones of 0. In figure [l] of this 
note like in figure 5 of |l[ the data can seem at first view well aligned on a straight line, 
i.e. decreasing with a simple exponential behavior. We will show now that this is just not 
true, and that these data clearly follow a non-exponential behavior. 

The fact that the decay is not purely exponential can be clearly seen, for example, by 
plotting only the points with r = 4 on the same semilogarithmic scale: we do that in figure 
|2[ Here it is clear that the points are not on a straight line. It is just the dispersion of the 
many points in figure [IJthat make them looking like being on a straight line: they are not, 
and figure |2| shows it clearly. 

We continue our analysis by trying to double check the conclusion of || leading to the 
functional dependence 



G(r,t) = 




We start by plotting in figure [5] the same data of figure versus ( |§y J 2 , together with the 
(very good) best fit to the form 
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Figure 3: G(r = 4,t) versus (|) 2 . The dashed curve is a straight line. 



G(r) exp (-g) . (7) 

(at fixed r as a function of t for the two parameters G and A). It is clear from figure |3] and 
from the quality of the best fit that this is the scale where the data follow a straight line, 
not the one of figure 

We repeat the procedure of figure |3] for different distances (finding always a very good 
best fit), and we fit the results to a power law (see figure |). We fit 

G(r) = r €l , A(r) = r t2 , (8) 

and for the best fit E\ = —0.54 and 62 = 1.47, in very good agreement with our initial 
Ansatz. 

Finally, in figure |5|, we show G as a function of r for our larger available time (t = 
819200), and the best fit to the form a exp (— b r% J . The best fit is again very good. It 
is interesting to note that, as a final confirmation of our point, the behavior of the function 
plotted in figure |5] is clearly not the one of a simple exponential, but that in a limited r 
region it can be mistakenly fitted with a simple exponential. 

We believe that these data and the data of [|IJ, when correctly analyzed, give strong 
evidence in favor of a Replica Symmetry Breaking like dynamical behavior of 3D spin 
glasses, and that they show beyond any doubt that the droplet model cannot not describe 
the numerical data one obtains for the q — q correlation function G(r,t). 
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Figure 4: G(r) and A(r) versus r and the best fit (dashed lines) to a power law (see 
text for more details). 
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Figure 5: G(r, t = 819200) versus r. The dashed line is our best fit to -2j-exp (—b ■ ri 
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